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Abstract 

Control of quantum dissipative systems can be challenging because control variables are typically part 
of the system Hamiltonian, which can only generate motion along unitary orbits of the system. To transit 
between orbits, one must harness the dissipation super-operator. To separate the inter-orbit dynamics 
from the Hamiltonian dynamics for a generic two-dimensional system, we project the Lindblad master 
equation onto the set of spectra of the density matrix, and we interpret the location along the orbit to 
be a new control variable. The resulting differential equation allows us to analyze the controllability 
of a general two-dimensional Lindblad system, particularly systems where the dissipative term has an 
anti-symmetric part. We extend this to categorize the possible purifiable systems in two dimensions. 

1 Introduction 

Recent decades have seen the application of mathematical control theory to quantum systems in both physics 
and chemistry, as technological advances have allowed for greater precision in manipulation of these systems 
[T] [2] [3] [4] [5j . One particular area of interest is the possible construction of quantum computers, which have 
the power to perform algorithms not accessible to conventional computers. A major experimental obstacle 
to any implementation of such a computer, however, is the decoherence of the system under influence of the 
environment. While much progress has been made on the control of closed quantum systems [B] [7] [S] , work 
on open quantum systems has proved more challenging [9] [10] [11] [12] [13] . 

One important issue is that controls are nearly always in the form of Hamiltonian operators. This type of 
control is unable to directly affect the purity of a state [T3] or transfer the state between unitary orbits. To 
control purity one must use the dissipative dynamics to move between orbits. To this end, we wish to derive 
a differential equation that captures only the inter-orbit dynamics, and collects the remaining dynamics 
(along the orbits) into a new control variable. This can be done if we assume arbitrary control over the 
Hamiltonian dynamics. The resulting differential equation can tell us how the location along the orbit affects 
the motion between orbits. In this paper, we show that this can be done for two-dimensional systems subject 
to Lindblad dissipation, and the formalism can be extended to consider the problem of purification (for 
related work in two dimensions, see [9 [10 ). The preliminaries of Lindblad dissipation are outlined in section 

2, and the projection of the Lindblad differential equation onto the set of unitary orbits is discussed in section 

3. In section 4, we analyze the controllability of this equation for various choices of system parameters — in 
particular, the case where the anti-symmetric part of the dissipation is non-zero. In section 5, we present a 
theorem that specifies necessary and sufficient conditions for purifiability. 
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tC.R. is at the Department of Physics, University of Windsor, ON, N9B 3P4. Canada, rangan@uwindsor.ca. Research is 
supported by NSERC, Canada. 



2 Preliminaries 



The state of a closed quantum system is described by a norm-one vector in a complex Hilbert space that 
evolves according to the Schrodinger equation: 

^(t)) = - iH \i,{t)) 

In order for the norm to be preserved, the Hamiltonian operator H must be Hermitian. An open quantum 
system, on the other hand, is described by a trace-one, positive-semidefinite operator p on the Hilbert space, 
known as the density operator (or density matrix when working in finite dimensions, as we shall). The 
interpretation of this matrix is the following: an eigenvalue of p is the probability that the system is in 
the corresponding eigenstate. Since the matrix is trace-one and positive-semidefinite, these eigenvalues are 
non- negative real numbers that sum to one. A state in the closed system becomes a rank-one projection 
operator \ip) (4>\^\ The Schrodinger equation, when extended to the density matrix, becomes the von Neumann 
equation: 

±p(t) = [-iH,p(t)] 

Certain relevant quantities are invariant under the von Neumann equation. The density matrix at any 
time can be written p(t) = f/(i)jo(0)i7 _1 (i), where U(t) is unitary. Since matrices at different times are 
similar, the eigenvalues are constant. The purity of the system, which is defined to be yjtr(p 2 ), is also 
invariant since it is the 2-norm of the vector of eigenvalues. This has implications for quantum control. Since 
control variables typically appear in the Hamiltonian only, the control dynamics cannot directly alter the 
probabilities, or purify the state (i.e. achieve a purity of one). 

However, a system that interacts with the environment will have non-Hamiltonian dynamics. In general, 
this will be an integro-differential equation, but if one assumes the dynamics depends only the present state 
and not its history (i.e. the Markovian condition) and there is not explicit time dependence, the resulting 
differential equation is the Lindblad equation|15)|16j: 

-\p(t)) = [-iH,p(t)]+J2[LjpL] 2 {L)L ]lP }\ (1) 

The Lindblad operators {Lj} can be taken to be traceless, as adding a multiple of the identity al to Lj is 
equivalent to adding an operator \(aLj — aLt) to the Hamiltonian. An alternate equation, known as the 
Lindblad-Kossakowski equation, chooses a basis {lj} of the set of traceless n-dimensional matrices that is 
orthonormal relative to the inner product (A,B) = tr(A^B): 

-\ P (t)) = [-iH,p(t)\ + ^ (hpA - Mh'P}) ( 2 ) 

j,k=l v J 

where the coefficients a^-fe's form a positive-semidefinite matrix, known as the Gorini-Kossakowski-Sudarshan 
matrix. 

A Lindblad operator can be thought as a stochastic jump with recoil. Under the influence of one Lindblad 
operator, a pure state IV'XV'I m time St becomes a mixture of two states, Mi\ip)(ip\M\ -\- M 2 \ip) (ip\M\ . Here, 
M\ = L\fSt and M% = I — LSt. In other words, lip) jumps to the state , 1 =L\ip) with probability 

y/ (4>\L't L\ip) 

^{ip\tfL\ip)5t. This is a jump because the St appears in the probability only, meaning the destination 
state does not approach the original state as St — > + . Conversely, the second state in the mixture is 
— ^ ^|^t_L|i/>)<5t ^ ^ J ^)\ % l J )i wri ich is infinitesimally close to the original state. In other words, depending 

on L, there may be an infinitesimal recoil needed to compensate for the jump process. When L^L is a 
multiple of the identity (for example, when the Lindblad operator is a multiple of a Pauli matrix), the second 
state reduces to the original state \ifj), so that the jump is recoil-less. 

lr The bra-ket notation prescribes that a vector be written as \a) and its dual as (a|. Inner products are written {a\b) and 
outer products (or rank-one matrices) |a}(6|. 
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3 Projection of dynamics in two dimensions 



As mentioned in the introduction, control of open quantum systems typically involves control variables in the 
Hamiltonian. Hamiltonian operators, however, can only move states along unitary orbits, and not between 
orbits. The goal of this paper is to isolate the between-orbit dynamics for a generic two-dimensional system 
under Lindblad dissipation. Our starting point is the following control system: 

j t P= \- iu i a iiP] + \ Yl a 3k\VjPVk-^{vkO-j,p}\ (3) 

j—x^y.z j ^k—x ,y ,z 

where {o~j : j = x, y, z} are the Pauli matrices. The controls {uj} are unbounded and may take any value in 
R. Note that we have chosen our set of control Hamiltonians to span su(2). In other words, we can make any 
unitary operator up to a non-physical phase difference, and therefore we can move between any two states on 
a given unitary orbit in arbitrary time. We are neglecting any drift Hamiltonian H Q = c I + Y]j— T y z CjO~j, 
since the component along the identity matrix does not contribute to the dynamics, and the components 
along the Pauli matrices can be treated by re-calibrating the control variables: Uj — > Uj — Cj. 

The density operator can be written in terms of the Pauli matrices: p = \{I + ^2j=x v z n j°~j): where the 
rij's are components of the Bloch vector, such that i.e. n%. + + rr z < 1. Substituting this expressions into 
the equation ([3]), we get: 



1 fajL a . = J^l-iujo-j, ^n k a k ] + i ^ a jk Wj,(Tk] + a jkM^j^k ~ \Wk<rj,cri}) 

j j,k jk jkl 

j,k ~ jkl 

The Pauli matrices obey the relations 



[o-j,a k ] = 2iy^ei<r l 
i 

{<7j,(7fc} = 2Sj k I 



OjO-iOk - -{(TkO-j, ai} = Skicrj + 5jiUk — 2S jk ai 



1 

2 

Using these relations, the Lindblad-Kossakowski equation above becomes 



l j,k,l 3^,1 jl 

If we define bi = J2^ a jk^jki, and a^ k — a i h + a "i ; we have 

l j-k^l I jl 



In vector notation, we can write: 

tin 



b + tixn+(A t > -tr(A b )I)n (4) 
at 



where A s is the matrix with elements af 4 . 



Now we want to decompose this equation into dynamics along and between unitary orbits, p has eigen- 

2 



values ^tF, where r := \n\ and eigenvectors 



l fa):= ^| 1) + _a |2) 

Note the spectra correspond one-to-one with the values of r, the Bloch radius. It follows that the unitary 
orbits are concentric spheres, except for the completely mixed state, which corresponds to the point r = 0. 
So we can parametrize the orbits by r, which lives on the closed interval [0, 1], and characterize the motion 
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along orbits with the unit vector n — n/r. We must be careful with respect to the innermost orbit however. 
n is not defined there, which means that the differential equations which we will derive for r and n will have 
solutions that exist for finite times (those solutions correspond to trajectories of p that pass through the 
completely mixed state). 

Since r 2 = n ■ n , 2r% = 2n ■ ff and therefore %=h-%. So: 

^ = h-b + n ■ (u x n) + h- (A s - tr{A s )I)n 

The middle term vanishes, the first term is constant in r and the third is linear in r. We can write: 

^=n-b + r(h-(A s n)-tr(A s )) (5) 

To find the ODE for n, we use n = rh, which gives ^ = — ^n). So we get: 

^ = 2u x n + - (b - (b ■ h)h) + (A s - n ■ (A s n))h (6) 
dt r 

Our goal here is to view equation ^ as a control ODE where h is the control. This view requires that we 
have full control over n, and we claim that we do, in terms specified by the following lemma. 

Lemma 3.1. Let S be the sphere centered at the origin with radius one, let B be the associated closed ball, 
and let B* be the closed ball with the origin removed. Let n(t) be a piecewise differentiable function from a 
time interval [0, T] onto S such that the corresponding solution r(t) of equation ^ is contained in the interval 
(0,1]. Then there are piecewise continuous control functions u x (t), u y (t) and u z (t) such that equation 
has the piecewise differentiable solution n(t) = r(t)h(t) on B* . 



Proof. First re- write equation (pi): 



u x h = 1 ( — — -(b—(b - h)h) - (A s - h ■ (A s n))h 
2 \dt r 

Any equation of the form x x a = b, where a • b = 0, has solution x — a x b. It follows that we can choose 
the controls to be: 

u(t) = h x - ( — - -(b— (b- h)h) - (A s - h ■ {A s h))h 
2 \dt r 

= \ \ x A — jtj™^ x ^~ x ^ A ' Sfl ^ 

Since n(t), n(t) and r(t) are piecewise continuous, so is u(t). □ 

Note that the prescription for u(t) is unbounded as r — > because of the middle term. This is because 
the system cannot approach the completely mixed state from any direction: when n = 0, t? is fixed to be b 
regardless of the controls u(t). 

We finish this section by writing down an alternate version of ^ in terms of the eigenvalues of A s , which 
allows us to specify a given system in terms of six real parameters. Let a± > > 03 be the eigenvalues of 
A s . Let {bj : j = 1, 2, 3} and {nj : j = 1, 2, 3} be the components of b and n relative to the intrinsic axes 
of A s (whereas the subscripts x, y and z denote the components relative to the eigenvectors of the Pauli 
matrices). This gives: 

,3 3 

The six parameters obey the following inequality, which arises from the positive semi-definiteness of A: 

a\b\ + aib\ + a 3 bl < 4aia 2 a 3 (8) 

The positive semi-definiteness of A also ensures the positive semi-definiteness of A s , so we also have ax, a 2 .a 3 > 
0. 
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4 Controllability analysis 

For a fixed r, the right-hand side of equation ^ can be seen as a map from S 2 , the set of available controls, 
to the set of possible values of r. Since this is a smooth map from a compact set to R, the image should be 
a closed finite interval. To analyze the controllability of we define functions /m(?~) and f m {r) to be the 
right and left endpoints, respectively, of this interval. That is, fu{r) is the maximum possible rate at which 
r can increase, and f(m) the minimum, for a given value of r. It is clear that (|7j) is controllable on a closed 
subinterval of (0, 1) if fu > and /,„ < everywhere on the subinterval. To steer between two points 
and r/, we choose our controls so that r(t) = fw (r(i)) if < r/, or r(t) = f m {r(t)) if > rf. 

Some properties of fu and f m can be gleaned from inspection of the differential equation, which we 
collect into a proposition: 

Proposition 4.1. If /m0") : = sup{f(r)} and f m (r) '■= inf{f(r)} 7 



I- Jm and f m are non-increasing. 

2. lim r ^ + fia(r) = \b\ and lim r _>.o+ fm(r) = -\b\. 

3. /m(1) < 0. 

4- fm(r) < for all r and system parameters. / ra (r) 
(which requires that all aj 's and bj 's are zero. 



= for r > only for the trivial where a\ = 



5.1fb has non-zero magnitude, /^(r) has an isolated intercept tt € (0, 1]. 

Proof. 1. If a control vector h* achieves the maximum r at r — r* , then choosing that control for all 
r < 7'* can only achieve a larger or equal f, since the coefficient of r in the differential equation, 
Y]j—l ^(l — rij), must be non-negative. Similarly, if a control n* achieves the minimum at r — r* , then 
choosing that control for all r > r* can only achieve a smaller or equal r. Furthermore, if a\ and a-i are 
positive, the coefficient of r cannot be made zero, so in this case, we can strengthen "non-increasing" 
to "decreasing". 

2. As r —> 0+, the linear term in ^ can be neglected, and we must extremize b ■ ft. The range of this is 
clearly 

3. Since r cannot exceed one, r| r =i < 0. 

4. Non-positivity follows from 1) and 2). If a\ > 0, f can be always made negative by choosing n = (0, 0, 1). 

5. Non-zero b implies that at ai and a 2 are positive, which means that /m is strictly decreasing on (0, 1). 
This, together with 2) and 3) imply the existence of r^- 

□ 

Corollary 4.2. Ifb is nonzero, there is an interval (0, rx), which we call a trap, inside of which the system 
is controllable. Outside of the trap, on [rr, 1], the system is one-way controllable; that is, can be steered 
to rf in finite time if and only if rf < n. 

Proof. The statements in the proposition imply that f m (r) < < /jj(r) on (0, r*r), in which case we can 
steer Ti to Tf > r, by choosing the control that satisfies f = /m(^/) provided rf < r^. Conversely, to steer 
r*i to Tf < r i: we can choose the control that satisfies r = f m ( r f )- On the interval [r T , 1], /^(r) < 0, so r.; 
cannot be steered to > r{, but can be steered to rf < Ti by choosing the control that satisfies f = f m (rf), 
which must be negative. □ 

In the case that |6| = 0, there is no trap: f < for all r, and in fact we can say that 

- r(ai + a 2 ) < r < -r(a 2 + a 3 ) (9) 

where we can achieve the upper and lower bounds by choosing n to be (±1, 0, 0) and (0, 0, ±1), respectively. 
In the case that a 2 = a 3 = 0, the decay of r may be halted, but otherwise r will decay exponentially to 
zero at a rate above or equal to a 2 + a 3 . It is evident, then, that the presence of an asymmetric part in the 
dissipative term (represented by b) significantly enhances the possibility of control. 
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In order to calculate /m and f m for given r, we can use the method of Lagrange multipliers. In some 
cases, we can solve the resulting equations analytically, but in general one must find the roots of a sixth-order 
polynomial, so we must resort to numerics. Before considering the general case, we will look at a particular 
case that can be treated analytically. We consider the possibility that a two-level system can undergo one of 
two processes represented by the raising and lowering operators <7+ and <7_ at rates a+ and a_, respectively. 
If one constructs the Lindblad equation using this scenario, and expresses it in the basis of the Pauli matrices, 
one finds that a\ = 0,2 = \ a+ ~ a ~ ][ , a 3 = 0, 61 = 62 — and 63 = a + — a_. The fact that b has only one 



non-zero component simplifies the equations so that we can treat the system analytically. 

If we apply the method of Lagrange multipliers to the right-hand side of |7| and set bi = 62 = and 
ai = a 2 , we get: 

2ra\ri\ = 2\ri\ 
2ra±n2 — 2Afi2 
63 = 2An 3 

71? + TU, + 7lo = 1 



where A is the Lagrange multiplier. This has the following solutions: 

n = (0,0, ±1) 



ni,n 2 , 



2a\r 



(10) 
(11) 



where n\ and rt 2 in (111 can be any pair that satsifies the normalization condition. Solutions (11) do not 
exist for all r, since the magnitude of 713 must not exceed one. They exist only on [J^rj, 1]- To determine 
which solutions correspond to /m and / m , we substitute back into (j7j). Solutions (10) give 



and solutions (11) give 



\h\ s 



Aair 



2a,\r 



(12) 



(13) 



We can easily conclude that f m (r) 



2a\r. Furthermore, the right-hand side of (13) is greater than 



or equal to those of (12), but since it has a limited interval of definition, we have: 



f M (r) 




2a ir , re(0,M) 
-rai, re 



(14) 



It happens that rj- in this case coincides with the point at which Jm switches between (10) and (11), i.e 



rx = j^j- This is not a general phenomenon, however: if 0,3 > 0, the switching point and the trap radius 
would not coincide. Fig. [I] depicts these solutions for a\ — (12 — 10 and 63 = 12. 

More generally, one can perform this analytical treatment in the following cases: (1) if 6 has one non-zero 
component, (2) if b has two non-zero components, and the corresponding a,j's are equal, and (3) b has three 
non-zero components, and a± — 0,2 = 03. If the system does not fall into any of those three categories, 
Lagrange multipliers lead to either a fourth-degree polynomial in A (technically solvable, but inordinately 
messy) or a sixth-degree polynomial (generally not solvable). The fourth-degree polynomial arises in the 
cases (1) b has two non-zero components but corresponding o^ 's are not equal and (2) b has three non-zero 
components and a\ — 0,2 > 03 or a\ > 0,2 = 03. The sixth-degree polynomial arises if b has three non-zero 
components and a\ > 0,2 > 03. 

In those cases, we can find the real roots of the polynomial numerically. Then we can compute the 
corresponding values of f, choose the maximum and minimum values, and assign the values to /m and f m . 
In fig. [2] /m and f m are shown for a particular system that required solving a sixth-order polynomial. We 
have computed the curves for 10,000 points apiece, tt can be found by numerically interpolating Jm- For 
the case depicted in fig. [2j tt was computed to be 0.544387876644064 (to machine precision). 
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Maximum achievable dr/dt 




0.2 0.4 0.6 0.8 1 

r 



Figure 1: Maximum and minimum achievable dr/dt vs. r for a case that can be solved analytically. System 
parameters: a\ = a 2 = 10, a 3 = 0, bi = b 2 = 0, b 3 = 12. Solid lines represent /m and f m . Blue and purple 
indicate solutions (10) and (11), respectively. Dotted lines indicate where these solutions do not coincide 
with f M . 




Figure 2: Maximum and minimum achievable dr/dt vs. r for a case that must be solved numerically. The 
trap radius is where the maximum achievable dr/dt passes from positive to negative. System parameters: 
ai = 10, a 2 = 5, a 3 = 0.3, h = 0.15VO, b 2 = 0.9, b 3 = 3V6 
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5 Purifiable systems 

An important goal in quantum control is purification: the process of steering a mixed state to a pure state, 
which can be characterized by a purity \/tr(p 2 ) = 1. Alternatively, a system is pure if the leading eigenvalue 
is one, with remaining eigenvalues being zero. In terms of the above analysis, we say a system is purifiable 
if and only if the trap radius rx — 1- In other words, the function /m 0") has an isolated intercept at r = 1. 
This section is devoted to proving a theorem that characterizes the possible purifiable systems. First, we will 
use the following lemma: 

Lemma 5.1. r = at r = 1 if and only the system is in a state that is an eigenvector of all contributing 
Lindblad operators. 

Proof. Because the Bloch radius can be written r = A + — A_, where A+ > A2 are the eigenvalues of p, we 
can write r — (ip + \p\ip + ) — (■)/>_ |/o|V-'-) ■ Differentiating this, we get an alternative expression for r: 

dr • • 

- = (ii>+\p\i> + ) - (4>-\p\il>-) + ^+\p\4>+) - ty-W-) + - 

= A + «v>'+|V+) + (V+IV'Y)) - A_«V>-|V-) + (<MV>-» + W+\p\i>+) - (tM/#-> 
= {i> + \pty + ) - (i>-\p\ii>-) 

where in the last step, the normalization of the vectors makes the quantities in parentheses vanish. Now, 
if the dissipation is characterized by a collection of Lindblad operators {L^j's, which are not necessarily 
orthogonal we can use ([T]) to specify p: 

- = (v>+|H*MIV>+> - (rP-\[-iH,p}\^) 



-{^\L jP Lt^_) + \{^\L]L lP \i>_) + \{^\ P L)L^) 



The Hamiltonian terms vanish since they are diagonal elements of a skew-symmetric matrix. We are interested 
in f when r — 1, so insert p — \ip+)(tp+\. We get: 

§ = £ (&+\L j \i>+){i>+\L]\i>+) ~ ^ + \L)L^ + ) - ^_|£ 3 # + >^+|it|^_) 

3 

If we insert the identity operator between and Lj in the middle term, we get the expression: 

J = -2^|(^|L,|^ + )| 2 
3 

For r to vanish, we need |( - 0_|L ) |?/'+)| 2 to vanish for each Lj. This is only possible however if \ip+) is an 
eigenvector of each Lj, since otherwise Lj\tp + ) would have some component in the \ip—} direction. This proves 
the lemma. □ 

This leads to the following theorem: 

Theorem 5.2. A two-level Lindblad system is purifiable if and only if one of the following characterizations 
hold: 

• There is one Lindblad operator, and it is singular. 

• There is one Linblad operator and it is non-singular with non- orthogonal eigenvectors. 

• There is no more than one singular Lindblad operator and any number of non-singular operators. All 
share a common eigenvector. 

• There are any number of non-singular Lindblad operators that share a common eigenvector. 



Proof. We are required to show two things to prove a system is purifiable: /m(1) = 0, and ai > 0. ai > 0. 
The latter ensures that Jm is strictly decreasing rather than constant in r. When combined with the former 
condition, this implies that /m is positive for all r < 1, and therefore controllable. 

It follows from the lemma that a system is purifiable only if all contributing Lindblad operators share a 
common eigenvector, or else /m(1) will be strictly negative. This is only a necessary condition however and 
not a sufficient one, since the condition implies only that /m(1) = 0. We also require that 02 > 0. So consider 
the case = 0. This implies that 03 and b are also zero (due to (JsJ) ) , so that A has only one non-zero entry 
in its natural basis. We claim that A in this form corresponds to a non-singular operator with orthogonal 
eigenvectors. It is a rank-one real positive matrix, and therefore can be written A = 53ij= x y z m i m j f° r some 
real 3-vector to. When one diagonalizes the Lindblad equation however, this results in a single Lindblad 
operator L = ^2j =x y z mjO~j- This operator is Hermitian and traceless, however, so neglecting the trivial 
zero operator, it is non-singular with orthogonal eigenvectors. 

In other words, as long as the system obeys the terms of the lemma, and does not consist of a single 
Hermitian operator, the system is purifiable. The first two cases in the theorem cover the remaining possible 
single-operator cases. The remaining two cases can be seen by noting that two singular operators cannot 
share eigenvectors, since they have only one (we consider two operators that are multiples of each other to 
be essentially one process). The third case covers the possibility of one singular operator: it has only one 
eigenvector, and that eigenvector must be shared with the other non-singular operator. The fourth case in 
the theorem covers the possibility of no singular operators but more than one non-singular operator. Note 
that the non-singular operators in the third and fourth cases need not have non-orthogonal eigenvectors. □ 

6 Conclusions 

We have shown that the inter-orbit dynamics of a controlled quantum system can be isolated from the intra- 
orbit dynamics by projecting onto the set of spectra of the density matrix. If one makes certain assumptions 
about the controllability of the system along the orbits, the position of the system along the orbit can be 
viewed as a new control variable, since the intra-orbit dynamics can be made arbitrarily faster than the inter- 
orbit dynamics. In two dimensions, we have derived a equation describing this inter-orbit dynamics, where 
the new control is the normalized Bloch vector, and the most general Lindblad system can be described 
by six real parameters: three describing the symmetric part of the dissipation, and three describing the 
anti-symmetric part. 

We have analyzed the controllability of a general two-dimensional system under Lindblad dissipation, 
particularly for dissipation with non-zero anti-symmetric part. For systems of this type, there exists a trap, 
or a subinterval of the state space where each state is reachable from any other, but from which states may 
not escape. The size of this trap can be calculated analytically for certain simple cases, but in general must 
be calculated numerically. We have shown how this can be done using the method of Lagrange multipliers, 
and shown results for a particular generic system. 

Furthermore, we have applied this formalism to categorize the set of purifiable systems. A necessary 
condition for purifiability is that all Lindblad operators share a common eigenvector. To strengthen this to 
a sufficient condition, one must eliminate the case of a single Hermitian Lindblad operator. 

The immediate direction of future work is to apply this formalism to three and higher dimensional systems. 
It is well-known that the structure of density matrices is richer and less well-understood than the case for 
two dimensions [T^j. For one, the set of pure states no longer constitutes the boundary of the set, but a 
(measure-zero) subset of the boundary. Currently we are studying how our formalism translates to higher 
dimensions and what obstructions are posed by the richer geometry. Furthermore, we would like to know 
whether the set of purifiable higher-dimensional systems can be categorized as it has been done in this paper. 
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